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I. INTRODUCTION 

Charge-carrier diffusion coefficients in semiconductors 
have been studied for decades, owing to their impor- 
tance as parameters in the operation of semiconductor 
deviceSfi"~ indeed, diffusion currents sometimes greatly 
exceed drift currents. The diffusion of minority carriers 
may differ from the same carriers' diffusion coefficient 
when they are majority carriers?^ and so must be deter- 
mined independently. The difference may be particularly 
important when the majority carriers are degenerate. 

Many devices contain degenerately-doped regions, 
such as the base of a bipolar junction transistor. Under 
conditions of high injection both majority and minority 
carriers may be degenerate, as occurs in the active layer 
of a diode laser. The diffusion coefficient of degenerate 
minority carriers, however, has not been explored. 

In this work we measure the diffusion of minority 
electrons by measuring Da, the coefficient of ambipolar 
diffusionii When an excitation such as the absorption of 
light creates equal numbers of excess electrons and holes, 
each species briefly diffuses according to its own density 
gradient. As the positive and negative charges become 
spatially separated, the resulting electric field, the Dem- 
ber field, causes the more mobile electrons to be held back 
by the less mobile holes. The resulting conjoined motion 
of the two species is ambipolar diffusion. Because the de- 
gree of charge separation is small, ambipolar diffusion is 



Electronic mail: lcweber@scu.edul 



usually described under an approximation of local charge 
neutrality. 

Due to the fiexibility of optical measurements, the 
use of photoexcited carriers has emerged as a prominent 
method for measuring ambipolar diffusion. ^^-'^^ However, 
the density of photoexcited carriers can strongly influence 
the measured diffusion coefficients)^ ' "'^^ Experiments on 
n-type semiconductors give a consistent result: Da is 
either an increasing or a flat function of photocarrier 
density. For instance, recent work on n-GaAs quantum 
wells has shown Da (or the related "ambipolar spin dif- 
fusion") to be increasing in lightly doped wells^^ and flat 
in nominally undoped wells.™*!- Bulk p-type semicon- 
ductors, however, have been less studied, and the results 
vary: Paget et alJ^ recently observed ambipolar diffusion 
in GaAs that strongly decreases with excitation density, 
while Zhac^'' saw no density dependence in silicon. Here 
we report measurements of heavily Zn-doped InP and 
flnd Da flrst to rapidly increase and then to level off. 

Electron-hole scattering has been shown to influence 
the dependence of Da on photocarrier concentration^i^ 
particularly in p-type quantum wells at low doping;^ 
However, for bulk material we flnd that even neglecting 
density-dependent scattering rates, we can harmonize all 
three behaviors seen in our and others' data — increase, 
decrease, and fiat — in a simple calculation following the 
approach of van Roosbroeck.1 We assume local charge 
neutrality, but account for the influence of carrier den- 
sity on the charge susceptibility. The regime of increasing 
Da that we observe is unique to heavily-doped samples 
at low temperature. In this regime Da nearly equals the 
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electrons' diffusion coefficient. 

Our ultrafast measurement allows us to separate sev- 
eral physical processes by their time-scales. The slower 
part of our signal reveals ambipolar diffusion. The faster 
part, lasting just a few picoseconds, reveals the effects of 
electron cooling and trapping. In this regime we are able 
to measure a very rapid ambipolar diffusion of holes and 
hot electrons. In keeping with our interpretation of the 
longer-time data, this rapid diffusion is due primarily to 
the decreased charge susceptibility of hot electrons. 




II. METHODS t[ps] 



To measure the diffusion and dynamics of photocarri- 
ers we use an ultrafast transient-grating method. A pair 
of "pump" laser pulses are simultaneously incident on the 
sample. The two pulses are non-collinear and interfere; 
their absorption excites photocarriers in a sinusoidal pat- 
tern with wavelength A and wavevector q — 27r/A. By 
locally modifying the index of refraction, the photocar- 
riers create a "grating" off of which time-delayed probe 
pulses diffract. As photocarriers recombine and diffuse, 
the grating amplitude decays at a rate 
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Here D is the diffusion coefficient, and tq is the lifetime of 
spatially uniform excitation. Measuring T{q) at several 
q determines D. We measure the diffracted probe am- 
plitude in a reflection geometry, improve the efficiency 
by heterodyne detectioufi^ and suppress noise by 95-Hz 
modulation of the grating phase and lock- in detection.^^ 

The pump and probe pulses come from a mode-locked 
Ti:Sapphire laser with wavelength near 800 nm and rep- 
etition rate of 80 MHz. The pulses are focused on the 
sample to a spot of 145 /im diameter. As the laser's flu- 
ence is varied, the probe pulses are always a factor of 6 
weaker than the pump pulses. At 800 nm InP has an 
absorption length of order 0.3 /im and reflectivity of 0.3 
(Ref. 113), so at our highest fluence each pair of pump 
pulses photoexcites electrons and holes at a mean den- 
sity of Hex ~ 2 X 10^'' cm~^. 

Implicit in Eq. [T]is the assumption that n, and there- 
fore Da, is independent of position. The pump pulses' in- 
terference excites densities from in the grating's troughs 
to 2nex at the peaks. However, the recombination time in 
^»-type InP is 33 ns (Ref. [l^, while our interpulse spac- 
ing is 12 ns, so a steady-state population of photocar- 
riers accumulates. Diffusion spreads these carriers uni- 
formly, raising the troughs to about 2ncx and the peaks 
to 4nex, and making Da roughly position-independent.^" 
The mean value 3nox corresponds, at our highest fluence, 
to an electronic quasi- Fermi energy Ecp/ks ~ 380 K. 
The steady-state photocarrier population assists in cool- 
ing hot electrons to the lattice temperature, since carrier- 
carrier thermalization is more rapid than thermalization 
through phonons. 



FIG. 1. Curves showing the time-dependence of the transient 
grating's amplitude for high-fluence (2.1 fij/crn^) excitation. 
Values oi q are 1.26 x 10* cm~^ (diamonds), 2.51 x 10'' cm~^ 
(squares), and 4.53 x 10* cm"'^ (circles). Lines are least- 
squares fits to Eq. [2] with B — 0. Curves decay more quickly 
at high q, indicating carrier diffusion of 34 cm^/s. 



The pump pulses create a carrier population within 
about 0.3 nm of the sample's surface. As the carriers 
diffuse inward some leave the probed region, causing the 
grating's diffraction efficiency to decay. This decay pro- 
cess may contribute to I/tq but will not influence the 
measurement of D because the inward diffusion is insen- 
sitive to the grating's in-plane wavevector q. 

All data we report are on a p-type sample)^ Zn- 
doped to a room-temperature carrier concentration of 
4.5 X 10^^ cm~'^, which corresponds to Ep/kB ~ 190 
K. It is oriented (100) and has room-temperature resis- 
tivity p — 2.65 X 10~^ cm and Hall mobility hh — 53 
cm^/Vs. Except where noted all transient- grating data 
are at 15 K. 



III. RESULTS AND DISCUSSION 

Figure [T] shows typical time-dependence of the sig- 
nal S{t) diffracted from the transient grating after high- 
ffuence excitation. S{t) rises slowly, reaching its full am- 
plitude only after several picoseconds, then decays ex- 
ponentially. The time evolution of the transient-grating 
amplitude consists of three distinct pieces, well described 
by the equation 

Sit) = -^e~*/^™ + Be~*/^p-k + Ce-*/^f»". (2) 

Here Tj-isc is the characteristic time of the curve's rising 
component, typically a few picoseconds; l/rfaii is the rate 
of the grating's relaxation. The term B exp (— t/rpcak) is 
zero except under very low-fluence excitation (Section 
HUB 2p , in which case Tpcak has a value intermediate be- 
tween Tj-isc and Tfaii- We begin in Section IIII Al by dis- 
cussing the dynamics of the falling component, and take 
up the rise and the peak in Section IIIIBI 
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FIG. 2. (a) Grating decay rate 1/r vs. at 15 K. Falling 
component: 1.2 /xj/cm^ (circles) and 0.15 /ij/cm^ (squares). 
Lines are fits to Eq. [1] The higher slope of the line at high flu- 
ence indicates significantly faster ambipolar diffusion than at 
low fluence. The rising component of the signal (diamonds, 
1.2 fiJ/cra^) diffuses still more quickly, (b) Diffusion as a 
function photoexcited density, as measured by excitation flu- 
ence. Diamonds: diffusion of the rising component at the two 
highest fluences. Circles: ambipolar diffusion of the falling 
component. Lines connect the points. At low fluence Da 
is electron-dominated and rises rapidly; at higher fluence an 
apparent saturation of Da is a crossover to hole-dominated 
diffusion. 



A. Ambipolar DifTusion 

By the time that the transient-grating signal begins 
to fall exponentially, the photoexcited electrons and 
holes have largely equilibrated to the lattice temperature. 
Electrons and holes are constrained to move together, re- 
sulting in ambipolar diffusion. 

The lower two curves in Fig. [UJa) show the decay rate 
l/Tfaii as a function of q^. The agreement with Eq. [T] 
reveals that carrier motion is diffusive. Strikingly, exci- 
tation at a fluence of 1.2 /uJ/cm^ results in electrons and 
holes that diffuse much more quickly than they do when 
excited at 0.15 ^J/crn^. Figure [5]Jb) shows the variation 
of Da', at low fluence it rises swiftly with the density of 
photoexcited carriers, before saturating at high density. 

We explain the observed variation of Da in terms of 
the electron and hole mobilities fie and fih, the back- 
ground hole density po, and the density n of photoexcited 
electrons and holes — which increases linearly with laser 
fluence. Under the assumption of local charge neutral- 
ity, one obtains the well-known result^ for the ambipolar 
diffusion coefhcient: 



Let /ichom be the electrons' chemical potential, and de- 
fine the electronic susceptibility as Xe = '^^/'^Mchom' 
fine hole susceptibility Xh analogously. Then diffusion is 
related to mobility byj^i^ 
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with an analogous equation for holes. These equations 
are just the Einstein relation, written in a form that is 
valid at all temperatures. In three dimensions Xe must be 
evaluated numerically^ but in the non-degenerate limit 
Eq. [5] reduces to the familiar expression Dg = /iefc^r/e. 
In the degenerate limit 
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where N is the density of states. 
Combining Eqs. [3] and HI 
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We can write the last factor as 1/x* = 1/xe + l/x/i- As n 
increases with fiuence, fie will decrease (as will fih, more 
slowly) ; but this decrease is typically quite weak,-^ If we 
treat fie^h as constants, Eq. [6]can explain our observation 
that Da rises and then saturates. In fact Da will vary 
non-monotonically, giving varied behaviors that match 
the variety of trends observed in Da vs. fluence liS^J^iiiii^ 

To gain insight into Eq. |6l we evaluate it first in sev- 
eral limiting cases, then numerically. When n <^po, at 
low fluence, Da = nfie/ex*- This result differs from the 
electrons' unipolar diffusion coefficient De by the factor 
Xe/x*j but at low fluence, Xe ~ X*- For sufficiently high 
n electronic degeneracy causes Da to increase with n, as 
observed. To reach high n while maintaining n <^ pq, as 
at our lowest fluences, requires a heavily-doped sample. 

At fluences much higher than those used in this ex- 
periment, n^po, the diffusion is proportional to the 
hole mobility and increases with n: Da — nfih/ex* ■ The 
crossover from electron-dominated to hole-dominated 
ambipolar diffusion results in the apparent saturation of 
Da seen in our data [Fig. HUb)] and calculations [Fig. 
[3fb)]. If the ratio h = fie/fJ-h is sufficiently large, the 
crossover region may be approximated by both n <C po 
and nfie :> PofJ-h- fn this case Da ^ pofih/ex* , giving 
Da OC n~^/^ at low temperature. In fact for InP b is of 
order 10 (Ref. IH), so the crossover is narrow enough to 
appear flat rather than decreasing. 



Da = 



{^) 1. Numerical evaluation of Eq. 



De and Dh are the electron and hole diffusion coefficients; 
tJe and tJ/i are the conductivities due to only the electrons 
or holes, respectively. 



To further elucidate our and others' results on p-type 
semiconductors, we evaluate Eq. [5] numerically. To avoid 
the complication of density-dependent mobilities, we plot 
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TABLE I. Parameters used in calculating Da/fJ-h, taken from Refs.~ '"'"°'"'"^'^^ . ml and are density-of-states effective masses 
in units of the free electron mass. Calculations of x assume parabolic, isotropic bands. 



the ratio Da/ fJ-h for a few representative values of b, using 
the InP electron and hole masses shown in Table ID The 
results for lightly-doped p-type material at 15 K and 300 
K appear in Fig. [3ja). For both temperatures, at low 
photocarrier density Da/ = bksT/e is a constant. It 
decreases in the crossover region, then increases in the 
hole-dominated, high-fluence regime. (By contrast, Eq. [B] 
shows that for n-type samples Da never decreases with n, 
consistent with measurements in nominally undope d^^i^^ 
and Si-doped'^- GaAs quantum wells.) 

Figure [21[b) shows Da/^h calculated for a sample 
doped as heavily as ours. For low temperature, even 
at low excitation density the degeneracy of photoexcited 
electrons causes Da/ fih to increase with n; this behavior 
does not occur at high temperature, or in lightly-doped 
samples at any temperature. The crossover regime man- 
ifests as the saturation seen in our highest-fluence data. 
Beyond this saturation, in the hole-dominated regime 
Da/ fJ-h continues to increase. 

For heavy doping and 300 K [Fig. [2Ib), dashed lines], 
the calculated Da/fih does not increase at low fluence. 
Low signal precludes measuring our sample at 300 K, but 
we have measured Da at 80 K [Fig. HJa)]. It confirms 
our expectation by behaving intermediately between the 
high- and low-temperature predictions, increasing more 
weakly than at 15 K. 

Turning to others' data. Fig. lUJb) shows the evaluation 
of Eq. ini under the conditions of Refs. [l3 andfl^ With 
increasing fluence, these papers found Da to be unchang- 
ing and strongly decreasing, respectively. Accounting for 
the experiments' doping levels and material parameters 
(Table m, the calculations reproduce both trends. The 
figure also shows calculated curves for the specific condi- 
tions of our own measurements, using our best estimate 
of the mobility ratio, b = 10. The 80 K prediction for Da 
does increase weakly with n, as observed. 



B. Short-time dynamics 
1. Rise 

The transient grating's behavior at short times, while 
electrons are still hot, further supports the analysis pre- 
sented above. Figure [5] shows the signal for high-fiuence 
excitation. At higher q the curve not only falls more 
quickly, but also rises more quickly. A plot of l/rriso vs. 



10 



=! 0.1 

CO 



J 00 



lOr 



(a) 





1 


.J^^!^..^,.***-. _ 


















(b) I 



















0.01 



Fluence [n/p 



100 



FIG. 3. (Dual-log scales; color online) Numerical calculations 
for Da/fih, for T = 15 K (solid) and 300 K (dashed). For each 
panel, curves are 6 = 3 (higher) and 30 (lower). Each curve 
is normalized to its low-fluence value, (a) Lightly doped, 
Po = 10^'' cm~^. (b) Po = 4.5 x 10^* cm~^, corresponding to 
the sample measured in this work. 
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FIG. 4. (a) Measured diffusion vs. fluence at 80 K. Da, as 
expected, increases more slowly than at 15 K. The diffusion of 
the rise remains much faster than that of the fall, (b) Da/i-th 
normalized to its low-fluence value, calculated using parame- 
ters from Table [H Under the conditions of: our experiment at 
15 K (1) and 80 K (2); Zhao at 90 K (3) and 300 K (4); and 
Paget et al. (5). Calculated trends agree with those observed. 
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FIG. 5. High-fluence (2.1 /ij/cm'^) transient-grating curves. 
Values ol q are 1.26 x lO"* cm"^ (diamonds), 3.77 x lO" cm"^ 
(squares), and 5.28 x 10* cm~^ (circles). Lines are least- 
squares fits to Eq. [2] with B = 0. Higher-g curves rise 
more quickly, indicating that carrier motion is diffusive even 
at short times. The diffusion coefficient of the rising com- 
ponent, 110 cm^/s, is much faster than the diffusion of the 
falling component. 

[upper curve of Fig. ^a.)] reveals diffusion that is much 
faster (3.2 times greater) than the ambipolar diffusion of 
the fall. 

The increased ambipolar diffusivity of hot carriers 
arises from two factors. First, in polar, cubic semiconduc- 
tors hot electrons' diffusivity has been observed to exceed 
that of cold electrons by factors up to fivei^ Under opti- 
cal excitation holes are heated much less than electrons, 
so we suppose that /ie increases while Hh remains un- 
changed; then Eq. |6] shows that electronic mobility can 
increase Da by a factor no greater than 1.6, because in the 
crossover regime at high fluence Da is strongly influenced 
by hole mobility. The remaining increase comes from Xe, 
which depends on the electrons' density and tempera- 
ture. A modest electronic heating of about 45 meV, for 
instance, would make the electrons non-degenerate and 
increase 1/x* by a factor of two. 

The electrons' initially high temperature may also give 
the rise its distinctive shape: a sign opposite to the rest 
of S{t), and a short lifetime. Our pump and probe pho- 
tons have the same energy, and the diffracted signal S{t) 
arises from transient, local changes to the complex in- 
dex of refraction evaluated at the probe's energy. The 
pump pulse initially excites electrons to states above the 
conduction-band minimum. As the electrons cool they 
move to lower energy states that are no longer degener- 
ate with the probe-photon energy. We speculate that this 
cooling may be accompanied by a shift in the primary 
origin of S{t) — for instance from phase-space filling to 
bandgap renormalization or free-carrier absorption. The 
latter two mechanisms are expected to yield a sign oppo- 
site to that of phase-space fillingj^ 

The rise's lifetime of a few picoseconds reflects the elec- 
trons' rapid thermalization. This thermalization slows at 
high fluence [Fig. El^a)], in keeping with many,—"— but 
not all,'^^ prior observations. The slowing is attributed 
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FIG. 6. (a) As a function of excitation fluence, lifetime of 
the rising component at 15 K (circles), rising component at 
80 K (squares), and peak at 15 K (diamonds). Both features 
decay more slowly at higher fluence. Lines connect the points, 
(b) Ratio at 15 K of the coefficients A/C (circles) and B/C 
(diamonds) in Eq. (2] At higher fluence the rise and peak 
features are both suppressed relative to the falling component. 
Ratios are the mean over all measured q. (c) Low-fluence 
transient-grating signal (suppressed zero) at 15 K. Squares: 
0.3 /iJ/cm^, q = 1.76 x 10'' cm"^ Diamonds: 0.08 /^J/cm^, 
q = 1.26 x 10* cm~*. In the lower-fluence data a prominent 
peak-like feature precedes the slow decay. At slightly higher 
fluence the peak begins to blend with the slow decay. Lines are 
least-squares fits to Eq. [2] with B 0. (d) 15 K transient- 
grating curves (log scale) at 0.08 /xJ/cm^. Diamonds: q = 
1.26 X 10* cm"^ Circles: q = 4.53 x 10* cm~\ Despite the 
complicating presence of the peak, the diffusion of the fall is 
clearly visible, though only 2 cm^/s. 



to the hot-phonon effect:— hot electrons couple primar- 
ily to LO phonons. The LO phonons couple weakly to 
acoustic phonons, giving heat only slowly to the lattice 
but frequently back to the electrons. This effect is partic- 
ularly strong in InP because of the wide energy gap be- 
tween acoustic and optical phonon branches. '^'^ The slow- 
ing may be compounded by increased screeningS^^^S, or by 
laser heating, as seen by comparing the rise lifetime at 15 
K with that at 80 K. The rise at 80 K [Figs. SJa), HJa)], 
though longer-lived, behaves the same as at low temper- 
ature: at high fluence it diffuses much more quickly than 
does the fall. As at low temperature, the increased diffu- 
sivity owes more to the hot electrons' susceptibility than 
to their scattering rate. 



2. Peak 

In addition to the change in diffusion rate Da with 
excitation fluence, the transient-grating signal's shape 
changes. Most markedly, at fluences of of 0.3 /zJ/cm^ and 
lower, S{t) acquires an additional component at short 
times, which appears as a "peak," shown in Figs, ^c) 
and[n{d). Fits to the form of Eq. [2] then require B ^ 0. 

As fluence increases, the size of the peak decreases rel- 
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ative to the overall signal size [Fig. [6jb)]. This is consis- 
tent with the peak arising from the trapping of a fixed 
number (and therefore a decreasing proportion) of car- 
riers. Because the peak appears at low fluence, where 
signal is small, we cannot determine whether it diffuses, 
but we plot its lifetime (the mean over all q measured) 
in Fig. ^sl). The peak's lifetime increases with increas- 
ing fluence. This slowing may arise from the electrons' 
slower thermalization: the cross-section for trapping of 
high-energy electrons is lower than that for low-energy 
oncs."^^ We believe the peak does not disappear abruptly 
at high fluence, but its reduced size and its slower decay 
renders it indistinguishable from the main decay process, 
Tfaii, due to ambipolar diffusion. 

Most important to the present work, the existence of 
the peak does not introduce any ambiguity into our low- 
fluence measurements of Da- The ambipolar diffusion 
takes place during the signal's slow fall. As Fig. IH^d) 
shows, after the peak has completely decayed the slow 
fall is the signal's only remaining component. 



IV. CONCLUSIONS 

We performed transient-grating measurements of am- 
bipolar diffusion in p-type InP for a wide range of pho- 
toexcited carrier densities. Ambipolar diffusion is com- 
monly described by two "rules of thumb" : that Da is con- 
trolled by the minority carriersr'"— or by the less-mobile 
carriersr^"— in a p-type sample both cannot simultane- 
ously be true. In fact when n <C po, at low fluence. Da 
is always nearly Dg. Our measurements in this regime 
revealed diffusion that increases strongly with n due to 
the electrons' degeneracy. 

As n increases our measurements reach a regime of 
intermediate density where Da is proportional to hole 
mobility but strongly enhanced by the factor Xh/Xe- 
Here the ambipolar diffusion levels off, but the several- 
picosecond rise of the transient grating's amplitude re- 
veals the distinct and much faster ambipolar diffusion of 
holes and hot electrons. The rise diffusion reaches 110 
cup /s, due largely to the effect of heating on the elec- 
trons' susceptibility Xe- 

We presented a simple calculation that assumes local 
charge neutrality and neglects the density-dependence of 
carrier mobilities. This calculation reproduces the salient 
features of our data — Da^s initial rise and subsequent lev- 
eling off — provided one accounts for the electrons' de- 
generacy through Xe- It also reproduces the Da seen by 
ZhaoiS- and Paget et alJ^ in more lightly-doped p-type 
samples, which are level and falling, respectively, over a 
broad range of excited densities. 

At yet higher electronic density Eq. [5] (Fig. [3]) pre- 
dicts that Da will increase with n for any temperature or 
doping level. At such high fluence Da may be dominated 
by the less-mobile carriers^^ but this regime in p-type 
samples has yet to be observed. 

Previous optical measurements have not seen diffusion 



of degenerate minority carriers, even though — as is evi- 
dent from Fig. HKb) — optical experiments routinely ex- 
cite electronic densities corresponding to Ecp/kB > 300 
K, and experiments with amplified lasers reach densities 
far higher. Rather, electronic degeneracy is usually ob- 
scured by the crossover to hole-dominated diffusion that 
begins when nfie is of order pofj-h- Diffusion coefficients 
under conditions of heavy doping and high injection are 
important to the operation of many semiconducting de- 
vices. Our transient-grating measurement highlights the 
important role of degeneracy in determining diffusivity, 
even for minority carriers. 
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